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Statistical description of hydrodynamic processes for ionic melts is proposed with tak- 
ing into account polarization effects caused by the deformation of external ionic shells. 
This description is carried out by means of the Zubarev nonequilibrium statistical 
operator method, appropriate for investigations of both strong and weak nonequi- 
librium processes. The nonequilibrium statistical operator and the generalized hy- 
drodynamic equations that take into account polarization processes are received for 
ionic-polarization model of ionic molten salts when the nonequilibrium averaged val- 
ues of densities of ions number, their momentum, dipole momentum and total energy 
are chosen for the reduced description parameters. A spectrum of collective excita- 
tions is investigated within the viscoelastic approximation for ion-polarization model 
of ionic melts. 
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I. INTRODUCTION 



Study of equilibrium and nonequilibrium properties of ionic melts remains actual from the 
viewpoint of experimentaU"— , tlieoreticaUi""— as well as computer simulation2ii22i"— investi- 
gations. They are very important because of a wide range of applications in chemical, metal- 
lurgical and nuclear technologies^^—. It is important to notice the papers^i^ where the binary 
distribution functions for a number of alkaline- chloride melts and some valence asymmetric 
ones were experimentally obtained for the first time by the method of neutron scattering in 
the systems with isotope replacement, developed by I.E. Enderby and colleagues. Dynamic 
structure factors of some ionic melts were received by means of an inelastic neutron scattering 
in worka^i"-. Moreover, diffusion, electroconductivity and viscosity coefficients were exam- 
ined experimentally^ as well. The recent X-rays inelastic scattering experiments on molten 
NaCl^, NaP and CsCf' have initiated study of collective excitations and dispersion laws 
of ionic melts. Theoretical investigations intended to explain observable phenomena in such 
systems were carried out on the basis of the kinetic equations for one-component plasma4^i^, 
the mode coupling theoryi^i^i^, the generalizedr^ and the extendeci^'^ hydrodynamics, the 
Zubarev nonequilibrium statistical operator (NSO) approacb^i^, the generalized collective 
excitations approacb2i>2^"-^>^ and others. It is important to note that in consequence of col- 
lective modes analysis for two-component plasma, five hydrodynamic modes (one heat, two 
sound and two mass and charge diffusion modes) and five relaxation modes were indicated 
based on the renormalized kinetics theory in papers^i^. Six of those modes describe prop- 
erties of a system as neutral and four modes as a charged one. Correlation functions and 
corresponding response functions of mass and charge densities, temperature and momentum 
divergency fiuctuations were considered in paper— using linear hydrodynamics equations, 
which takes into account thermoelectric and electrostriction effects. Statistical hydrodynam- 
ics of ionic systems were constructed in paper— with help of the nonequilibrium statistical 
operator method^i^ based on the extended set of the reduced description parameters, in- 
cluding microscopic densities of particles number, their momentum, total energy and densi- 
ties of generalized viscous stress tensor and energy flow. The obtained equations of extended 
hydrodynamics are valid for both strong and weak nonequilibrium processes. Moreover, in 
paper— this approach was reformulated on the basis of the generalized Fokker-Planck equa- 
tion for the collective variables functional, considering nonlinear hydrodynamic fluctuations 
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in ionic systems. The approach presented in paper— permitted to discover mutual influence 
of heat- viscosity processes in the time correlation functions "mass-mass", "mass-charge", 
"charge-charge" and their flows for weakly nonequilibrium processes in NaCl ionic melt. 
Theoretical study of dynamic structure factors, time correlation functions of longitudinal 
and transverse currents of mass and charge densities as well as transport coefficientsii""— 
obtained a qualitative agreement with the results of molecular dynamics calculation32Z"i2^ 
actively carried out in the mid 70s of the last century, starting from the paper—. In that 
time, the important results in investigations of structure and dynamical properties of ionic 
systems by means of the molecular dynamics method were obtained. In particular, the equi- 
librium binary distribution functions of ionic melts, calculated with the molecular dynamics 
method^S have a good agreement with the experimental data^'^. By means of the molecular 
dynamics simulations the dynamic structure factors for the model ionic melts of NaClM^ 
and RbBr^, time correlation functions of ionic velocities, diffusion and electroconductivity 
coefficients for ionic melts of NaCl, Lil, RbClM'^ and others were received. An interesting 
fact is that the spectrum of charge density fluctuations of ionic melts in the long-wavelength 
limit has a characteristic peak like a spectrum of longitudinal optic phonons of ionic crystals. 
This was shown using the molecular dynamics simulations in paper™ and afterward in the 
theoretical studiesi^""— . Moreover, it was confirmed by Ab initio MD calculations in conjunc- 
tion with the generalized collective mode approacb^i and in real experiments^. Experimental 
investigations^"-, MD computer simulations^i^^ and Ab initio MD simulations^ii^Si"^ reveal 
an importance of taking into account polarization processes caused by electronic structure 
of ions. On the other hand, the theoretical study of collective modes dispersion for NaCl in 
the generalized collective mode approach^ showed that using the hard spheres ionic model 
yields too high values of optic modes frequencies. It is worth to mention that in theoretical 
studies and MD computer simulations of the equilibrium and nonequilibrium properties of 
ionic melts an effective ionic Haggings-Majer— and Tosi-Fumi^^i^ potentials which do not 
take into account polarization effects were used. However, in the real melts outer electrons 
shells can be polarized and so must contribute in an effective interaction. Today polarization 
effects are described within the three models: polarization point dipole models {PPDMs), 
shell models (SMg) and fluctuating charge models (FCMg)^. In PPDMg an inducted point 
dipole is added to atomic or ionic positions, whereas in other two models dipole has a flnite 
length. In SMg dipole is presented via a pair of point charges, namely with a positive nucleus 
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and a negative shell, linked with a harmonic connection; whereas in FCMg point charges 
are fixed in certain positions and their values can fiuctuate. In particular, infiuence of po- 
larization effects on equilibrium properties of ionic melts Agl were considered in study^. 
In paper>2^ polarization effects are described within PPDMg model when studying dynam- 
ical properties by means of MD simulations. In order to take into account electron shells 
deformations in works^^i^ a shell model was proposed in which electron shell is presented as 
an external electron cloud that interacts with nucleus through the harmonic potential and 
repulsive potential at small distance^^. Based on this model, the Car-Parrinello Ab initio 
simulations in which the dynamics of ionic subsystem was reduced to the pseudo-dynamics 
of electrons wave functions within the density functional formalism were carried out. This 
permits observation of polarization effects caused by deformations of external electron shells 
at the Ab initio level. 

In present paper a statistical description of hydrodynamic processes in ionic melts with 
taking into account polarization effects caused by deformations of external electrons shells 
of ions is proposed. It is implemented by means of the Zubarev nonequilibrium statistical 
operator method^'^ that permits to study both weak and strong nonequilibrium processes. 
In the second section the nonequilibrium statistical operator and the generalized hydrody- 
namic equations taking into account polarization effects are obtained for the ion-polarization 
model of ionic melts, when the nonequilibrium averaged values of densities of ions number 
n"(r), their momentum p"(r), angular momentum s'^(r), total energy ^(r) and dipole mo- 
ment d"(r) are chosen for the reduced description parameters. The generalized molecular 
hydrodynamics equations for ionic melts in the case of weak nonequilibrium are obtained 
in the third section. Section 4 contains analytic expressions for the collective excitations 
spectrum in the limit k — )■ 0, w — on the basis of the viscoelastic model for ionic melts 
with taking into account polarization processes. 

II. HAMILTONIAN OF THE SYSTEM. NONEQUILIBRIUM 
STATISTICAL OPERATOR 

Let us consider an ion-polarization model of ionic melts that classically describes ionic 
subsystem with taking into account polarization effects. We assume that both positively and 
negatively charged ions can be polarized, though in case of Nal ionic melt, only negatively 
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charged jodine ions are polarized. In polarization processes related with electron transition 
between orbitals of an atom that turns to a positive ion and a negative atom of melt that 
turns to a negative ion, dipole moments of ions are induced. In consequence of interactions 
like "induced dipole - induced dipole" , "induced dipole - ion" and motion dynamics it leads 
to induced momenta of dipoles that in turn cause their rotational motion. Therefore, ions 
with induced dipole moments besides translational degrees of freedom possess rotational 
ones due to polarization processes. All these degrees of freedom of ions in melts must be 
accounted in the full Hamiltonian of the system. The latter can be presented in the following 
form: 

a j=l V 'I / 

which includes kinetic (translational and rotational energies) and potential parts; pj denotes 
a momentum vector, rria is a mass of ion of species a, wj stands for an angular velocity and 

^ y 

Ja means an inertia tensor of the j-th polarized ion (ionic dipole) determined relatively its 
center of mass. Uion we introduce as follows^^: 

U.on = 2 E E '^'^^(Ir^^D -EE - 2 E E (2) 

a,b i,j=l a j=l b j=l 

+ E E ^ + E E f-birtj)^rijd„ 

a j=l a,b i,j=l 

where indicates that i ^ j when a = b. $a{,(|ry|) is an ion- ion interaction potential: 

$a.(|rd) = ^f^, (3) 
I ^tj I 

Za, Zjj are the valences of the ions of a corresponding species, e denotes an electron charge, 
(Ta, <Jb are the radii of the ions of a corresponding species, dj means a dipole momentum of 
an ion with a distorted external shell 

b i=l ' 

aj is a polarizability of the j-th ion in an electrical field 

E, = + EJ, (5) 

b t^j 

Xdi-Tij) 1 



b i^j 



\>-tJ\ I i I 



\rij\ = \ri — Tjl is a distance between ions. Damping dispersion functions /abdrjjl) can be 
determined according to paper—. It is important to notice that interaction of polarized 
ions has a central-asymmetric character. Obviously, effects related with induced rotational 
degrees of freedom are expected to be small comparing to the induced dipoles processes in 
the general hydrodynamics picture of ionic melts. 

Nonequilibrium states of the ion-polarization model of ionic melts are described by the 
nonequilibrium statistical operator p{x^;t), which satisfies the Liouville equation 

^p(x^;t) + zLjvp(x^;t) = 0, (6) 

where iL^v stands for the Liouville operator, corresponding to Hamiltonian ([T]): 



= E E — ■ 7^ + i^A) ■■^] (7) 



( Q d-d d 



dj = 1^ is unit vector describing spatial orientations of ionic dipole. 

To solve the Liouville equation ([6]) we use the Zubarev nonequilibrium statistical operator 
method, in which solutions are searched according to the N. Bogolyubov's idea of a reduced 
description of nonequilibrium processes based on the set of observed variables (P„(r))*. Using 
this method, the solution of equation (jH]) can be present in a general form with taking into 
account the projection: 

p(x^;t) = PreAx'^-^t) - I e'^''-'^TUt',t){l - Prei{t'))iL^prei{x'';t')dt', (8) 

J — oo 

where 

Trel{t',t) = exp+ |- (1 - Prelit"))lLNdt"'^ 

is the evolution operator with taking into account Kawasaki-Gunton projection Prei{t). Pro- 
jection operator depends on a structure of relevant statistical operator prei{x^;t) 

PreiW = (prM'';t) - [ ^r^^^^ (P„(r)) ^ Sp(p') (9) 

V / dr^P4^Sm.ir)p') 
d{P^{r)Y 



and possess the following properties 

Prel{t)Prel{t) = Prel{t), (1 - Prelit)) Prelit') = 0, 

Prel{t)p{t) = Prel{t), Prel{t) Prel{t') = Prel{t)- 

Relevant statistical operator Preiit) is received from the conditions of an informational en- 
tropy extremum at fixed values of the reduced description parameters (P„(r))* including the 
normalization condition Spp^eK^^j^) — 1- Within Gibbs approach one can obtain 



(10) 



Preiix^'-t) = exp -Y.J drFn{r;t)Pn{r)^ , 

=lnSp expj-^y drF„(r;t)P„(r)|. 

is the Massieu-Planck functional, F„(r;t) are the Lagrange multipliers, which are de- 
termined from the self-consistency conditions 

(Pn(r))* = (Pn(r))*e. (11) 

and thermodynamic relations 

Here, S{t) denotes entropy of nonequilibrium state of the system determined according to 
Gibbs 

sit)^m + T.^n{t){Pn{r)y. (13) 

n 

We will consider hydrodynamic state of ionic melt within the formulated ion-polarization 
model. For its description the averaged values of densities of ionic number fi"(r) = 
(5(r — Tj), their momentum p"(r) = Ylf=iPj^i^ ~ ^j)' angular momentum s"(r) = 
Ej=i^j^ir - r,), total energy e{r) = EaEf=i (it + 1^ ' w,w.)5(r - r,) + 
Xlj^i ^ion(^i)'^('^ ~ ^j) along with an induced dipole moment d"(r) = X^^i dj5(r — r^), 
which are observable variables and satisfy corresponding conservational laws, can be chosen 
for the reduced description parameters. For such a set of variables prei{x^] t) can be written 
down as follows: 

Preiix^'-t) = expj - $(t) - J drP{r;t) [s{r) - ^ v"(r;t) • p"(r) (14) 
- J2 ^'^(r; t)n''{r) - E(r; t) ■ d'^(r) - ^ w'^(r; t) ■ s%r)] |, 

n n n ^ 



where the Lagrange multiphers (3{r; t), v"(r; t), z/°(r; t), E(r; t), w°(r; t) are determined from 
the self-consistency conditions 



= m)u (is) 

(p'^(r))* = {r{r))U 
{n^{r)y = {n\v))U 
{d\v)Y = mr))U,, 
(s'^(r))* = (s"(r))*,, 

and from the nonequihbrium thermodynamic relations ( IT2l) . ( fT3l) . They have the following 
meaning: /3(r; t) = l/kBT{r; t) is the inverse local temperature; v'^(r; t) is the mean value of 
hydrodynamic velocity of ions; w"(r;t) denotes mean value of angular velocity of polarized 
ion, z/'^(r;t) = /i«(r;t) + Ilh^q^ + ^-w-(r^)w^(r;t) , ^ ^ Z,e<^(r;t) stands 

for the electrochemical potential; /i"(r; t) is the chemical potential of ions; ip{r; t) denotes 
scalar potential of electric field E(r; t) induced by ions and dipoles of the system. Electric 
field E(r; t) = (E(r))* satisfies the averaged Maxwell equations: 

Vx(E(r))* = -l|(B(r))*, (16) 
Vx(H(r))* = i|(D(r))* (17) 



+ Tfe^(P"Wr + 5:-ld".V(p"(r))A 

\ a a / 



V ■ (B(r))* = 0, (18) 



V • (D(r))* = A7rlj2 ^ae(n"(r))* + J]] d" • V(n"(r))* j 

\ a a / 



(19) 



Where, microscopic electric E(r) and magnetic H(r) fields and corresponding inductions 
D(r), B(r) satisfy the microscopic Lorenz-Maxwell equations. Known integral relations 
between (D(r))* and (E(r))* as well as between (B(r))* and (H(r))* determine spatially 
inhomogeneous dielectric function e(r, r'; t, t') and magnetization x(r, r'; t, t') which describe 
polarization processes in the system. Acting by the operators {1 — Prei{t)) and zLjv on Prei{t) 
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in (IHl) we obtain 



(1 - Prei{t))iLr,Prei{t) = { " / drf3{r; t)/,(r; t) (20) 
+ drf3ir;t)E{r;t)r,{r;t)-J2 j t^r/3(r; t)w'^(r; ^/^(r; t)|p.e«W, 



where 



/,(r;t) = (l-P(t))2Ljv£(r), (21) 

i;{v-t) = {l-P{t))iLMv\^), 

4"(r;t) = (l-P(t))^L;vn"(r), 

/,"(r;t) = (l-P(t))^L;vd"(r), 

j:(r;t) = (l-P(t))^LArs'^(r), 

are the generahzed flows, and I^{r]t) = 0. P{t) denotes the generahzed Mori projection 
operator which in this case has the following structure 

Pm = {A)l, + [ rfr||i^(l(r) - (^(r))*) (22) 



V^(d'*(r))* ^5(s'*(r))* J 

and properties P(t)(l - P(t)) = 0, P(t)Pn(r) = P„(r). 

Taking into account (1201) we obtain the nonequilibrium statistical operator of the ion- 
polarization model of ionic melt 

p{t) = Preiit) + j' e<'~'">Trei{t,t') j rfr |/3(r, t')4(r; t') (23) 

- 5^ /3(r; t') v"(r; t')/;(r; t') - /3(r; t')E(r; tO/,"(r; t') 

a a 

- ^ /3(r; t')w%r; t')i:{r; t') \prei{t')dt\ 



in which the generalized flows of energy density Ie{r;t), momentum density Ip{T;t), angu- 
lar momentum density /"(r;t) and dipole moments density I^ir^t) describe the dissipative 



processes in the system. Using the NSO (l23l) we can obtain generahzed hydrodynamics equa- 
tions for the reduced description parameters (P(r))* = {(n"(r))*, (p"(r))*, (s"(r))*, (d"(r))*, 
{i{r)y} within the ion-polarization model of ionic melt. We present them in matrix form 

d_ 

di' 



;(P(r))* = (P(r))t,, + I dr' I e^(*-*')^,,(r, r'; t, t')F{v'- t')dt\ 



(24) 



where P(r) is the column vector, P(r) = zLjvP(r), F{r';t') = {— /3(r'; t')z/"(r'; t'), 
-/3(r';t')v"(r';t'), -/3(r'; t')w'^(r'; t'), -P{r';t')E{r';t'), -f3{r';t')} is the column vector of 
the nonequilibrium thermodynamic parameters, ipii{r,r';t,t') is the matrix of the general- 
ized transport kernels (memory functions): 

""O 



rel 



'f>pp ipps <ppd <Ppe 

<f>sp 'f'ss <f>sd <Pse 

<^dp 'pds 'f'dd <^d£ 

(fep (pes Ved Vee 



(25) 



iPpp{v,v']t,t') 



"ftp vtp 



Vpp ^pp 



(26) 



{r,r';t,t') 

^;l{r,r'-t,t') = {i;{r;t)Trei{t,tVy-t'))i^ (27) 

are the generahzed transport kernels describing viscous processes, herewith, and ippp~ 
define generalized coefficients of viscosity caused by translational motion of positively and 
negatively charged polarized ions (a, b = {+, — }). 



V9,^(r,r';t,)f:') 



J (r,r';i,i') 

^ti^y-At') = (/r(r;t)T.,,(t,t')/.^(r';t'))L 



(28) 



(29) 



are generalized transport kernels describing viscous ionic processes caused by rotational 
motion of polarized ions, y?^^"*" and ipj~ define the generalized coefficients of angular viscosity 
of positively and negatively charged polarized ions. 



ipddir,r';t,t') 



"fid "Pm 



"Pdd "Pdd 



(30) 



J (r,r';t,4') 
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(^^^(r, r'; t, t') = (/,«(r; t)TUt. t'))U (31) 

denote the generalized transport kernels that describe transport processes of dipole moments 
of polarized ions of corresponding species, where 

12{v-t) = (1 - P(t))zL;vd"(r) = (1 - P{t)) f -__n^(r) + J] (a;, • d,)5(r - r,) j(32) 

and n^(r) = E74d,P/(r-r,). 

(p,,{v,v'-t,t') = {I,{r;t)Trei{t,t')Ie{r'-,t'))ii, (33) 

is the generalized transport kernel of total energy, which determines the generalized coef- 
ficient of heat conductivity of ionic melt within the ion-polarization model. The matrix 
elements (ppp, (ppd, (ppd, <^pe, 'f'de in (l25|) describe cross dissipative correlations between mo- 
menta, dipole moments and total energy flows of the system. The obtained nonequilibrium 
statistical operator fl2^ and the generalized hydrodynamics equations ^IM together with 
the set of the generalized Maxwell equations f ll6p - (IT^ are valid for both weak and strong 
nonequilibrium processes in ionic melts with taking into account polarization processes. Set 
of equations fl2^ is unclosed and describe viscous, heat and polarization processes. The 
system of hydrodynamic equations can be significantly simplified and becomes closed for a 
weakly nonequilibrium processes, when the thermodynamic parameters -F(r; t) and the re- 
duced description parameters {P{r)y slowly vary in time and space and slightly deviate from 
their equilibrium values -Fo(r), (P(r))o respectively. In following section we will consider a 
weakly nonequilibrium processes in ionic melts within the ion-polarization model. 



III. GENERALIZED MOLECULAR HYDRODYNAMICS EQUATIONS 
FOR IONIC MELTS WITH TAKING INTO ACCOUNT POLARIZATION 
EFFECTS 

In the case when the nonequilibrium thermodynamic parameters F{r; t) slightly deviate 
from their equilibrium values -Fo(r), the relevant statistical operator ( !T4|) can be expanded 
in deviations 6F{r;t) = F{r]t) — -Fo(r) with taking into account linear terms only. Then, 
excluding F(r; t) from the relevant statistical operator by means of the self-consistency 
conditions (ITSl) one can obtain 

pUt) = Po ( 1 + Y.SPik;t)^-\k)P{k)] . (34) 
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Here, po is the equilibrium statistical operator of ionic melt, SP(k; t) = (-P(k; t))*— (P(k; 0))o, 
P(k) is a column vector whose elements are the Fourier-components of the reduced descrip- 
tion parameters {P(k) = / e^^''P{r)dr} = {n"(k), p"(k), s"(k),d''(k), ^(k)}. ^-\k) is the 
inverse of the matrix of equilibrium correlation functions $(k) 



i.(k) 



Here, 



$nn(k) 











$nd 


$ne 






















l> 
















^dd 


^ds 










^ed 






nn 


(k) 


$+- 
nn 


(k)" 




$-4 
nn 


(k) 


$~ 

nn 


(k) 



(35) 



(k) 



(36) 







is the matrix of static structure factors of ionic subsystem $^„(k) = (fi"(k)n (— k)) 
Sabi^), where (. • .)o = / dV^ ■ ■ ■ po- 

with ^;;{k) = (p'^(k)p»(-k))o, and 

*s.(k) = 



$--(k) 



(37) 



<J^i+(k) 
$.T(k) 

whose elements ^^"(k) = (s"(k)s"(— k))o. Similarly, 



(38) 



^dd(k) 



^^dH^) <^>.V(k) 
*.7(k) <f,T(k) 



(39) 



where $J^^(k) = (d"(k)d^(— k))o are the equilibrium correlation functions of Fourier- 
components of dipole moments densities for polarized ions a and b. 



$,,(k) = (e(k)f(-k))c 



(40) 



is the equilibrium correlation Kubo-like function of Fourier- components of total energy den- 
sity of ionic melt. Other matrix elements in (l35l) describe static correlations of variables 
n"(k), p"(k), d'^(k), £(k). 
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In the approximation (IMj) the nonequihbrium statistical operator has the following struc- 
ture 



where 



ab 



ab 



/,«(k) = (l-Po)^^ivd"(k), 



(41) 



/;(k) = (l-Po)^W(k), (42) 
I,{k) = {1 - Po)iLNe{k) 



are the generalized flows in a weakly nonequilibrium case. Pq denotes Mori projection 
operator 



PoA = (i)o + J2{A pW(-k))o<l'-^(k)P(k) 



(43) 



that possesses the following properties: Po(l — Po) = 0, PoP(k) = P(k). The set of the 
generalized hydrodynamics equations for ionic melt within the ion-polarization model with 
approximation fHTl) is now closed and can be presented in the matrix form 

|(P(k))*-zf](k)(P(k))*+ r e<'-''^^{k;t,t'){P{k)y'dt' = 0. (44) 



Here, 



iQ{k) = (tLNpik) p(+)(-k))o<i'"'(k) 



is the frequency matrix whose elements 



in(k) 



i^jip t^ns 



iVtsn iVtsd i^se 

iQdp iflds 





iQsrp iQss 



(k) 



(45) 



(46) 
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are the normalized static correlation functions. 



^{k;t,t') = {i{k)To{t,t')n^\-k))o^-\k) 





ippp (pps <^pd <^pe 

(fsp <^ss <^sd 'Pse 

<^dp 'f'ds <^dd <Pde 

. <^£p ifes <^£d <^e£ _ 



(47) 



(k;i,t') 



is the matrix of transport kernels (memory function), which describe a weakly nonequilib- 
rium transport processes including viscous, polarization and heat processes in ionic melts 
within the ion-polarization model, /(k) = (/^(k), J" (k), J^(k), /^(k)) is the column vector, 
/'^+^(— k) = (/^(— k), Jp (— k), /^(— k), /£(— k)) is the row vector. As can be shown, in the 
framework of the NSO method^i^, the time correlation functions of the basic set of the 
dynamic variables 



$(k;t) = (P(k;t)pW(-k))c 



(48) 



satisfy the equation (jUj) as well 

d ^ 



— $(k;t) -i(](k)$(k;t) + / e''^'-'>^(k-t,t')<l>(k;t')dt' = 



(49) 



which takes into account memory effects. In the Markovian approximation when these 
effects are negligible, in the limit k — 0, w — )■ and without taking into consideration cross 
dissipative correlations and rotational degrees of freedom, the set of equations fH9l) can be 
presented as follows: 



d ~ 



dt 



$(k;t)+T(k)<l>(k;t) = 0, 



where 



T(k) = -^fi(k) + <^(k). 



(50) 



(51) 



is the generalized hydrodynamic matrix. Here, 



iVt{k) = ik 



Unp 

Upn Upd UJpe 

Codp 

Usp 



(52) 



14 



is the corresponding frequency matrix and 

















4>dd 

0^ 



'ee 



(53) 



the matrix of memory functions, in which the wave-vector dependence is extracted. The 
system of equations fH^ permits to study the time correlation functions of partial dynamics 
in ionic melts. Though if the elements of the frequency matrix fj46|) can be calculated 
via equilibrium characteristics, in particular, the partial structure factors, the calculation 
of the transport kernels fH7|) is a well-known problem and it can be implemented only 
approximately with taking into account special feature of dissipative processes in the system. 
This issue accompanies the investigation of collective excitation spectrum for ionic melts with 
taking into account polarization effects, in particular in the hydrodynamics limit k — )■ 0, a; — )■ 
0. In the following section we use obtained above the generalized molecular hydrodynamics 
equations for ionic melts, in particular, within viscoelastic approximation, to investigate the 
spectrum of collective excitations in the hydrodynamic limit. 

IV. COLLECTIVE MODES OF IONIC MELTS. VISCOELASTIC 
APPROXIMATION. 

In order to study collective excitations within the ion-polarization model of ionic melts we 
will use viscoelastic approximation. We will not take into consideration total energy density 
as well as rotational degrees of freedom of induced dipoles and pass from the description 
in terms of partial dynamic variables n°'(k), p"(k), d'^(k) to the description based on total 
densities of ions mass n(k) = m+n"'"(k) + ^^-^^(k), total ions momentum p(k) = p'''(k) + 
p^(k), total charge ^(k) = Z+eh'^(k) + Z_eh^(k) and density of total current of charge 
j(k) = ^^p+(k) + ^l^p^(k). In addition we consider only negatively charged atoms to 
be polarized, i.e. negative ions possess induced dipole moments. Thus, we include density 
of dipole moment of negatively charged ions d^(k) to the set of the reduced description 
parameters. However, it is more convenient to use the variable ci(k) orthogonal to the 
variables ?7,(k) and ^(k): 



d(k) = d-(k) - A(k)n(k) - B{k)Q{k). 



(54) 
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Here, A(k) and -B(k) are determined by the relations 

A(k) = (d-(k)n(-k))o[<l>-^(k)]„„, + (d-(k)^(-k))o[<l>-^(k)] 



en 1 



B{k) = (d-(k)n(-k))o[$-^(k)]„,+ (d-(k)^(-k))o[$-^(k)],,. 

In this case, the generahzed hydrodynamic matrix has the following form 





(55) 
(56) 















k (j)pp 


-ikupi. 


kVpj 











—ikugj 




\<i^(t)jp 


—ik.ujjQ 


kVji 







3d 



\ 



(57) 



k^^dp - ^kwdp k^^dj - ikwdi k^^dd 
with the corresponding elements of the frequency matrix and the memory functions 
matrix constructed on the generalized flows /a(k) = (1 — Po)^-^Af^^(k), where a(k) = 
{n(k),p(k),^(k),j(k),d(k)}. 

0pp — is a normalized viscosity coefficient of ionic melt, 
0dd = 4>d-d- = D stands for normalized dipole diffusion coefficient of negatively 
charged ions, (pjj = a means an electroconductivity coefficient of ionic melt, where 



0J„ 



4tt 



is a squared plasma frequency of ions of species a. In 



(rp-)c 



the hydrodynamic matrix fl57j) the elements uis is calculated applying the |k| — )■ limit 
in a;/s(k). Taking into account (p^p^)o = hm 
lim|k|^o(p"(k)p"(-k))o = (n+ = ^, n 

of concentrations of negative and positive ions), we find for (pp)o the following result 



kKo(p+(k)p+(-k))o = \y y /o 

_ = ^ are the equilibrium average value 



lim|k|^o(p(k)p(-k))o 



■), or 



(pp)o = ^5 ''^ = n+m+ + n^m^ is 



-pp Nff /u ^ : 

the equilibrium value of the total mass density of ions of the system. Similarly, one can 
find $jj = (jj)o = lim|k|-s>o(j(k)j(— k))o = -^^p- Cross-correlation function is equal to 
$pj = (pj)o = lim|k|-i-o(p(k)j(— k))o = ^{Z^erij^ + Z^enJ), and taking into account the 
electroneutrality condition (Z+er;,+ + Z^enJ) = one can get that ^pj = $jp = 0. Con- 
sidering the structure of the normalized static correlation functions Uis and values of the 
averages (p^p^)o5 (p P )o) we find for uis the following results: 

m ~ 



np 



^nj — 0, 



pn 



UJ, 



P 
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QUI 



(58) 
(59) 



(^gp — 0, 



pg 



ngi 



(60) 



JQ 



471(3 



(61) 



gp 



due to condition of electroneutrality of the system. In the expressions ([58 



(EH) the quantities [^-^n, [^'%g, [^'^]gn, gg are the elements of the matrix $ 

inverse to the matrix of static correlation functions 















\ 


















































(62) 



\ <^dd / 

Whereas, = lim[k|^o '^'nn(k) = lim|k|^o{"^+5++(k) + m+m_[S+-(k) + 5-+(k)] + 

mlS-^{k)} = m^5++(0) + m+m_[S+-(0) + ^-+(0)] + mlS~'{0), where we use 5++(k), 
^^^(k), S'^^(k), S* (k) for the partial static structure factors. In particular, S'"*'+(0) = 
1 + n+/rfr(^^+(r) - 1) = ^xj, ^""(0) = 1 + n_ / dr(^^"(r) - 1) = ^x?, where ^^+(r), 
5'^^(r), xj, x^ are the equilibrium pair distribution functions of ions and their isothermal 
susceptibilities. Similarly, one can calculate "charge-density" $p„ and "charge-charge" 

Z+e[m+5++(0) + m_5+-(0)] + Z_e[m+5~+(0) + 
(Z+e)25++(0) + Z+Z_e2[5+-(0) + S-+m + (Z_e)25-(0). 



$pg correlation functions. We find ^ng 



Then the same way one can obtain the elements Upd, uj^p, oudj, oujd of the frequency matrix 



pd 



7" 



m 



d--A(0)-)[<I> 



\ddi 



(-d--^(o)-) |*-V 



'id 



Z_en. 
Z_en. 



■d- - 5(0) 



47r/3 / 



(63) 



(64) 



(65) 



(66) 



where [$ ^j^d is the element of the matrix $ ^ inverse to $, and "I'dd = lini|ki-j.o(d(k)d(— k))o. 
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In order to obtain the collective excitation spectrum of ionic melt with taking into account 
polarization properties we will consider the elements uis of the hydrodynamic matrix fl57|) 
and do not take into account dynamical cross correlations between flows of momentum, 
electric current and diffusive dipole flow which are described by the memory functions (ppd, 
(pdp, (pjd, <Pdj- Thereby, the hydrodynamic collective modes can be found as the solution of 
the equation 

^ z{k) —ikconp 

z{k) 

—ikcojn k'^cpjp —ikojjg z{k) + k'^cpj 



det 



—ikugj 





-ikupd 


-ikujd 



\ 



\ 



0. 



—ikudp —ikudj z{k) + k'^D J 



We will find the solutions as series in wave vector {k = |k|) 

z{k) = Zik + Z2k'^. 

Thus, one can get the set of equations for evaluation of the coefficients zi and Z2. After 
calculations we find that collective excitations spectrum contains two acoustic modes 



Zs± — i.iCsk V gk ] 



(67) 



two charge propagating modes 



Zch± — ^iCchk ^chk ; 



(68) 



and one diffusive polarization mode, related with dipole diffusion of negatively charged ions 

(^2(^2 _ \j \p 



Zd = -D- 



G 



en ne ^2 



Here, we use the definitions 



+ V(C4 - AG) 



1/2 



(69) 



(70) 



Cch — 



V2 



- ^{C^ - AG) 



1/2 



(71) 



that are sound and charge waves velocity respectively. = Cf + C2 + G^ + Gj, where 



^ ' ' ' ' '^'^ ' ' ' ' '^'^ ^dp^pd and C| — uJdjOOjd- G — + GfGl + C|C| 



Gi — UJnp(^pni G2 



^ngHn ~ KgKd^ whcrc wc introduce A;^fe = UacUJcb, for example A^^ 



^np^pQ- 



18 



In equations (l67I) - (l68l) Ts and Tch define the damping coefficients of corresponding waves 
determined by tfie relations 

1 



1 
2 



+ 



+ D-)- 



2(C2C2 - 2G) 



'pp 



Cl [Cl + Cl) - G 



+ 



'33 



Cl {cl + cl) - c 



D- 



Cl [cl + cl) -C- ClCl + K^k 



3P' ^Pj 



+ <f>P3 K + ^1p) + K K + Kj) 



(72) 



cn 2 



PP 



+ 



'33 



+ D-) 



1 



2(C2,C2 - 2G) 



'PP 



Clh {cl + cl) - G 



+ 



'33 



Clh {cl + cl) - G 



D 



Cl, {Cl + Cl) -G- ClCl + A^A. 



+ <PP3 (At + Al) + (A,', + K3) 



(73) 



Perturbation theory 

In the case when the cross correlations are sufficiently small one can get the simplest 
analytic expressions for collective modes. In order to do this we use the perturbation theory 
approach for collective modes developed in papers^>^. Its idea is to present the generalized 
hydrodynamic matrix as the sum 



T{k)=To{k) + 5T{k), 



(74) 



where the ffist term corresponds to the "bare" modes that can be easily obtained analytically. 
The second is considered as perturbation and its contributions into collective modes can be 
found from the expansion of the perturbation theory 



The following block matrix corresponds to the "bare" modes of the system 



To{k) 





-ikoj^ 






ikujnp 

pn k (f)pp 


















—ikuj 



S3 



-ikojjg k'^'Pjj 



\ 





k^D— J 



(75) 
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Solving the eigenvalue problem in "zero" approximation, we obtain the collective modes 
spectrum together with the corresponding eigenvectors: two conjugated acoustic modes 



z 



two charge propagating modes 



.(0) 



.2r(0) 



2;,^ = -^(±1,1,0,0,0), 



Xeh± = ^(0,0, ±1,1,0), 



and the diffusive polarization mode, related with dipole diffusion of negative ions 



.(0) 



0;^ = (0,0,0,0,1). 



Here, the corresponding damping coefficients are determined as follows: 



r 



(0) 
ch 



(76) 



(77) 



(78) 



(79) 



5T{k) 



(80) 



The cross correlations define the perturbation matrix 

/ 





-ikujjn k'^(f)jp 

—ikudp —ikudjO 

Now we can evaluate corrections for collective modes caused by the weak cross correlation. 
According to^"^, such corrections are equal to zero in the first order in perturbation, while 
in the second order they are determined by the expression 



\ 



\ 

-ikujpd 


-ikujd 

/ 



(0) _ ^(0) 



(81) 

where /3 runs over all possible values, except of a subclass. STai3 are the perturbation 
matrix elements in the representation of matrix To(A;) eigenvectors. STa/3 = {xa^Txji) = 
J2i j ^i.a^'^ij^j,l3- After calculation of the necessary elements 6Tai3 using the eigenvectors 
of matrix To(A;) and taking into account the results of "zero" approximation fl76l) - fl78|) we 
receive the final expressions for the collective modes. For two acoustic modes we obtain 



-Si 



^(0) I J2) 

^s± ^ ^s± 



±ikCi 1 



(82) 



CI 

2Cf 



(r 



(0) 



2Cf 
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for two charge propagating modes 



^ch± - + ^fh± (83) 

2 



±.,c.|i-H^ -.Mr™ 




and for diffusive polarization mode related with the dipole diffusion of negatively charged 
ions 

One can see, that the acoustic and charge propagating modes calculated within the per- 
turbation theory do not contain their cross-correlation contributions. This fact is caused by 
the asymmetry of the frequency matrix iVt due to the condition of electroneutrality of the 
system. Thus, the perturbation theory gives us the compact approximate expressions for 
the collective modes in which only the cross correlations related with the dipole moment are 
taken into account. On the other hand, when in expressions (1671) -( 168|) . that are explicit in 
— )■ 0, w — limit we neglect the correlations with the dynamical variable of the dipole mo- 
ment and the correlations between viscous and charge variables are also considered as small 
ones, we reproduce the results obtained within "zero" approximation of the perturbation 
theory (1761) -(177|). 



V. CONCLUSIONS 



We presented the statistical description of hydrodynamic processes of ionic melts with 
taking into account polarization effects caused by deformation of external electron shells of 
ions. It is implemented by means of the Zubarev nonequilibrium statistical operator method 
that enables to study both weak and strong nonequilibrium processes. As the result, the 
nonequilibrium statistical operator and the generalized hydrodynamics equations with tak- 
ing into account polarization effects are received within the ion-polarization model of ionic 
melts, when the observable values such as the nonequilibrium averaged values of densities 
of ionic number n'^(r), their momentum p"(r), angular momentum s''(r), total energy e(r) 
along with dipole moment d"(r) are chosen for the reduced description parameters. The 
generalized molecular hydrodynamics equations for ionic melts with taking into account po- 
larization effects in weakly nonequilibrium case are obtained as well. Based on them within 
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the viscoelastic approximation for ionic melts in the hmit /c — )■ 0, w — > for the spectrum 
of collective excitations we found two conjugated acoustic modes, two charge propagating 
modes and one diffusive mode, related with dipole diffusion of negatively charged ions. These 
modes where calculated using expansion in wave vector as well as perturbation theory. Here- 
with, in the limit of infinitesimal cross correlations from fl^7|) - fl68p the expressions of "zero" 
approximation can be reproduced. It is worth to mention that taking into consideration 
cross correlations described by memory functions (ppd, 4>dp, 4>jd and (pdj in both cases will 
renormalize the damping coefficients in propagating excitations and will not effect the shape 
of diffusive polarization mode. However perturbation theory yield approximate expressions 
for collective modes, it can be very useful as well. In particular, it permits to study the time 
correlation functions built on the basic dynamic variables fHS]) and investigate contributions 
caused by polarization effects in it within the generalized collective modes approach^. Per- 
turbation theory allow us to easily extend the model to thermoviscoelastic by taking into 
account the dynamic variable of total energy density in the region of small viscothermal 
correlations. 

Investigations of transport processes in multi-component melts LiF, BeF2, ThF^, UF4 for 
nuclear reactors with the uranium cycle'^° and the melts NaF, BeF2, LiF, ZrF^, Lz2-BeF4 
for the uranium and plutonium fuels^ are of great interest within the presented approach. 
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